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ABSTRACT 

The Fermi Gamma-ray Space Telescope (FGST) has opened a new high-energy window in the 
study of Gamma-Ray Bursts (GRBs). Here we present a thorough analysis of GRB 080825C, which 
triggered the Fermi Gamma-ray Burst Monitor (GBM), and was the first firm detection of a GRB by 
the Fermi Large Area Telescope (LAT). We discuss the LAT event selections, background estimation, 
significance calculations, and localization for Fermi GRBs in general and GRB 080825C in particular. 
We show the results of temporal and time-resolved spectral analysis of the GBM and LAT data. We 
also present some theoretical interpretation of GRB 080825C observations as well as some common 
features observed in other LAT GRBs. 
Subject headings: gamma-rays: bursts 
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1. INTRODUCTION 

Gamma-Ray Bursts (GRBs) originate from the most 
luminous explosions in the universe a nd more than 35 
years after their discovery in 1967 (Klebcsad el et all 
1973), many questions remain to be answered about 
their possible progenitors, the composition of the ultra- 
relativistic outflows that power them, and the dominant 
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emission mechanism for their prompt gamma rays. The 
Burst And Transient Source Experiment (BATSE) on- 
board the Compton Gamma- Ray Observatory (CGRO; 
1991-2000) made significant advances to the field, thor- 
oughly exploring the 25 keV - 2 MeV energy range with 
detailed population studies of the prompt gamma-ray 
emission. Burst spectr a were found to b e well described 
by the Band function (|Band et al.lll993ft . which consists 
of two smoothly connected power laws. It was under- 
stood, however, that observations of GRBs at higher 
energies were of crucial importance to resolve some of 
the open issues: constrain the bulk Lorentz factor of 
the outflow and the distance from the central source 
to the gamma-ray emission region, distinguish between 
hadronic and leptonic origins of the gamma-ray emission, 
and probe for signatures of Ultra High Energy Cosmic 
Rays (UHECR s) which cou ld be accelerated within GRB 
jets (see lBand et al.l [20091 . for a review of the prospects 
for GRB science with Fermi LAT). 

Constraints on the origin of the high-energy emission 
from GRBs are quite limited due to both the small 
number of bursts with firm high-energy detection and 
the small number of events that were detected in such 
cases. High-energy emission from GRBs was first ob- 
served by the Energetic Gamma-Ray Experiment Tele- 
scope (EGRET, covering the energy range from 30 MeV 
to 30 GeV) onboard CGRO. Emission above 100 MeV 
was detected in five cases: GRBs 9 10503, 910601, 930131, 
940217 and 940301 (|Dingusll2003l ). One of these sources, 
GRB 930131, had high-energy emission that was consis- 
tent with an extrapolation from its spec trum obtained 
with BATSE between 25 keV - 4 MeV (jSommer et al.l 
119941 ). In contrast, evidence for an additional high- 
energy component up to 200 MeV with a different tem- 
poral behavior to the low-energy component was discov- 
ered in GRB 941017 (in EGRET's calorimeter TASC; 
iGonzalez et"aH 120031 ). The high -energy emission for the 
latter GRB lasted more than 200 seconds with a sin- 
gle spectral component being ruled out. A unique as- 
pect of the high-energy emission in GRB 940217 was 
its duration, which lasted up to ~90 minutes after the 
BATSE GRB trigger, inc luding an 18 GeV photon at 
~75 minutes post-trigger (|Hurlev et al.| [l~994). More re- 
cently, the GRID instrument onboard Astro-rivelatore 
Gamma a Immagini LEggero (AGILE) detected 10 high- 
energy events with energies up to 300 MeV from GRB 
080514B, in coincidence wit h its lower energy e mission, 
with a significance of 3.0 a (Giulian i et al.ll2008f ). 

The Fermi Gamma-ray Space Telescope was launched 
on June 11 2008 and provides an unprecedented energy 
coverage and sensitivity for the study of high-energy 
emission in GRBs. It is composed of tw o instruments: 



the G amma-ray Burst Monitor (GBM; Meegan et al 



2009) and the Large Area Telescope (XAT: IAtwood et al 



20091) . The GBM covers the entire unocculted sky with 
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12 sodium iodide (Nal) detectors with different orienta- 
tion placed around the spacecraft and covering an energy 
range from 8 keV to 1 MeV, and two bismuth germanate 
(BGO) scintillators placed on opposite sides of the space- 
craft with energy coverage from 200 keV to 40 MeV. The 
LAT is a pair conversion telescope made up of 4 x 4 arrays 
of silicon strip trackers and cesium iodide calorimeter 
modules covered by a segmented anti-coincidence detec- 
tor designed to efficiently reject charged particles. The 
energy coverage of the LAT instrument ranges from 20 
MeV to more than 300 GeV with a field-of-view (FoV) 
of ^2.4 steradians. Note that the LAT effective area is 
still non-zero even as far out as 70 degrees off- axis which 
allows the detection of bursts with such high incident 
angles. As of June 1 2009, 9 GRBs have been detected 
by the LAT at ener gies above 100 M eV: GRB 080825C 
(IBouvier et al.ll200l . GRB 080916 C (lAbdo et alj|2 009a; 
iTaiima et aLll2008ft. GRB 081024 B (IQmodei et al.l l2008ft. 
GRB 081215A dMcEnerv et alj 12008ft. GRB~09021 7 
(jOhno et al. 2009a ), GRB 090323 jOhno et al.ll20 09b). 
GRB 090328 dMcEnerv et all 120091: iCutini et al.|] 2009). 
GRB 090510 (|Ohno et al.U2009d : IQmodei et al.ll2009i Tln 
this paper, we report the observations and analyses of 
gamma-ray emission from GRB 080825C, the first GRB 
detected by both the GBM and the LAT instruments. 
Section [2] will present the GBM and LAT observations 
along with the various methods used for data analysis, 
section [3] provides the results of detailed time-resolved 
spectroscopy, and section 0] discusses the theoretical in- 
terpretation of our observations and compares the prop- 
erties of this event to the ones observed in some other 
LAT GRBs. 

2. BURST DETECTION AND LOCALIZATION 

2.1. GBM observations 

At 14:13:48 UTC on August 25 2008 (T ), 
GRB 080825C triggered the GBM flight s oftware (trig- 
ger 241366429, [^mderHoriFeFS] [2001). On-ground 
analysis of the GBM data localized the burst at Right 
Ascension (RA, J2000) = 232.2°, Declination (Dec, 
J2000) = -4.9°, with a statistical uncertainty of 1.5° 
at the 1-cr confidence level. The GBM on-ground 
localization placed this GRB at ~60° from the LAT 
boresight at the time of the trigger, at the edge of the 
LAT FoV, where the effective area is a factor of ~3 less 
than on axis. 

The top two panels of Figure [1] show the background 
subtracted light curves (see section ET2"|) of the two bright- 
est Nal detectors (9 & 10) and of the two BGO detectors. 
The GRB exhibits a multiple peak structure with the 
two brightest peaks seen right after onset. The Tgo and 
T50 durations of the event (time duri ng which 90% and 
50% of the event flux was collected, cf. iKouveliotou et al.1 
119931 ) were estimated to be (8-1000 keV) ~27 s and 
~13 s, respectively. Emission in the Nal and BGO detec- 
tors becomes extremely weak after ~Tb + 25 s. However, 
emission is detected in the Nal scintillators up to ~35 s 
after the trigger time with a 3.6 a significance in bin (e). 

2.2. LAT observations 

The LAT events detected close to the GBM position 
around the trigger time are shown in the bottom panel of 
Figure[T](see details below). The LAT data show a count 



rate increase that is spatially and temporally correlated 
with the GBM emission. We have performed a detailed 
analysis of the significance of this detection. Details on 
this computation as well as of the LAT data selection, 
background estimation and localization are given below. 

2.2.1. Event selection 

Most of the events detected by the LAT instrument are 
cosmic rays that need to be distinguished from any source 
of 7-ray signal. The broad range of LAT observations 
and analyses, from GRBs to extended diffuse radiation, 
leads to different optimizations of the event selections 
which have different rates of residual backgrounds (mis- 
classified cosmic rays). The LAT background rejection 
analysis has been constructed to allow a nalysis classes 
to be optimized for specific science topics (jAtwood et al.1 
2009). In the case of GRB observations, the relatively 
small region of the sky as well as the very short time 
window allow the background rejection cuts to be re- 
laxed relative to an analysis of a diffuse source covering 
a large portion of the sky over longer periods of time. 
Indeed, the so-called 'diffuse' event class is most suited 
for studying faint sources (like diffuse gamma-ray emis- 
sion) with minimum background contamination. In the 
case of GRB 080825C, it was used to search for possi- 
ble high-energy afterglow emission up to 13 ks after the 
burst trigger (see section [3.1.2p . On shorter timescales, 
a significant improvement in signal-to-noise ratio can be 
obtained by increasing the effective area while keeping 
the background rate at a reasonable level. The so-called 
'transient' event class was developed for this specific pur- 
pose and is used for burst detection and localization. 

Using Monte-Carlo simulations as well as real data 
input, the event selection used for spectral analy- 
sis has been optim ized compared to previous analyses 
(jAbdo et al.l l2009aft that uses the 'transient' class for 
the prompt emission analysis. Indeed, we found that 
cuts can be loosened even further for typical burst dura- 
tion since background contamination is less of an issue 
for short time windows. For this purpose, a more re- 
laxed event class, so-called 'S3', has been developed to 
improve the LAT effective area at the expense of an in- 
crease in background rate. Pure effective area increase 
is on the order of ~30 % at 100 MeV and ~10% at 1 
GeV. For a typical GRB spectrum, the gamma efficiency 
above 80 MeV increases by a factor of 20% when using 
the 'S3' class (with a small dependence on spectral in- 
dex and incidence angle) when the all-sky-background 
rate increases from ~4.2 Hz ('transient') to ^5.2 Hz 
('S3'). Our study showed that 'S3' brings improvement 
over the 'transient' class in terms of signal-to- noise ratio 
(Signal/ y/ Signal + Background) above 80 MeV. This 
improvement depends on the brightness of the burst and 
is of the order of 15% for bursts with GRB 080825C 
characteristics. As a consequence, the 'S3' event class 
was used for the spectral analysis of GRB 080825C. Note 
that the lower energy threshold does not add any addi- 
tional systematics compared to the previously used 100 
MeV threshold (see section 1572")) . 

Finally, the event selection makes use of the spatial 
information around the best LAT localization. The 
LAT point-spread-function (PSF) has a strong depen- 
dence with energy as well as with the conversion point 
in the tracker. LAT events are thus separated into 
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FRONT (conversion in the upper part of the tracker) 
and B ACK (conversion in t he lower part of the tracker) 
events (|Atwood et al. 2009) for which separate response 
functions are provided. The region of interest (ROI) con- 
sidered in this analysis is energy dependent and based on 
the 95% containment radius (PSF95) and the 95% LAT 
error localization (Err95): 

ROI(E) = ^/PSF95(E) 2 + Err95 2 

To avoid large background contamination, a maximum 
size is set at 10 and 12 deg for FRONT and BACK events 
respectively. In the particular case of GRB 080825C: 
ROI(E < 200 MeV) is set to this maximum size and 
ROI{~ 500 MeV) = 2.9 and 4.0 degrees for FRONT and 
BACK events respectively. 

All events resulting from this selection process are 
shown in the bottom panel of Figure [T] 

2.2.2. Background estimation 

Since the number of events detected by the LAT is only 
15 (of which 13 have an energy >80 MeV and can can be 
used for spectral analysis), we need to carefully estimate 
the expected number of background events in order to 
calculate the statistical significance of our measurements. 
The background in the LAT data used for this analysis is 
dominated by cosmic rays (CRs), with a small contribu- 
tion from extragalactic and Galactic diffuse gamma rays. 
Because the earth's limb was ~85 deg from the GRB lo- 
cation, the contamination from earth-albedo events was 
negligible and consisted of events with very poorly recon- 
structed directions. The background rate is a function 
of many parameters and can vary by more than a factor 
of two, depending on the observational conditions. For 
example, there is a strong dependence of the CR back- 
ground rate on the geomagnetic coordinates at the loca- 
tion of the spacecraft. Furthermore, the background rate 
also depends on the burst position in instrument coordi- 
nates, because the LAT's acceptance varies strongly with 
the inclination angle. For these reasons, it is not straight- 
forward to estimate the expected amount of background 
during the GRB emission using off-source regions around 
the trigger time, since the spacecraft will have moved to 
regions of different geomagnetic coordinates, and the in- 
clination angle of the region of interest will have changed 
significantly All these effects are properly taken into ac- 
count to estimate the background rate for our specific 
observational conditions. 

Because the two background components (CRs and 
gamma rays) have different properties, they are esti- 
mated separately using two different methods. The 
amount of gamma-ray background from some direction 
in the celestial sphere depends only on the accumulated 
exposure in that particular direction. Therefore, this 
component can be estimated by simply scaling the num- 
ber of gamma rays detected in six months of LAT data, 
produced during normal science operations, by the ratio 
of the exposure of the GRB observation over the expo- 
sure of the six-month data set. Similarly to the above, 
the amount of CR background from some direction in 
the celestial sphere depends on the exposure in that di- 
rection. However, unlike the above, the CR background 
also depends on the geomagnetic coordinates at the loca- 
tion of the spacecraft at each instant of the observation. 



Because of the latter dependence, the CR background 
cannot be calculated the same way as the gamma-ray 
background. Instead, a Monte Carlo simulation of the 
GRB observation is performed, in which parameteriza- 
tions of the dependence of the CR background rate on 
the geomagnetic coordinates and on the inclination angle 
are used to estimate it for each second of the observation. 
These two methods will be described in detail below. 

The CR-background estimation was based on prop- 
erties of the LAT data extracted from a subset of six 
months of data selecting when the Galactic plane was 
far from the center of the LAT's FoV {\B\ > 70°). Dur- 
ing such observations, and for the transient and S3 data 
classes, the gamma-ray contribution to the detected sig- 
nal is negligible, and the detected events can be approx- 
imated as being cosmic rays. Specifically, for the tran- 
sient class the gamma-ray contamination in this subset 
of the data is about 10%, which consists of compara- 
ble amounts of galactic and extragalactic diffuse emis- 
sions and a negligible contribution from resolved point 
sources. For the S3 class, the gamma-ray contamination 
is comparable, although slightly lower. After extracting 
the dependence of the CR-background rate on the geo- 
magnetic coordinates at the position of the spacecraft, 
we used this dependence to calculate the all-sky CR- 
background rate for each second of the GRB observation. 
Then, a corresponding number of simulated events were 
generated with directions in instrument-centered coordi- 
nates, off-axis 9 and azimuthal angle <p, selected to match 
the observed distribution in the same subset of the six 
months of LAT data as used above. These coordinates 
were then converted to equatorial coordinates, using the 
instrument's pointing information and were added to a 
skymap. The simulation procedure described above was 
repeated hundreds of millions of times. The individual 
skymaps generated at each one of the iterations were then 
averaged to obtain a single skymap that showed the ex- 
pected amount of CR background from each direction of 
the celestial sphere for the specific GRB observation. 

As mentioned above, the gamma-ray component of the 
background was estimated by rescaling the number of 
gamma rays detected during six months of LAT data 
(without applying a cut on the galactic latitude of the 
LAT pointing direction this time). First, a CR back- 
ground map corresponding to the six-month period was 
created by following the procedure mentioned in the pre- 
vious paragraph. Then, an all-particle (both CRs and 
gamma rays) map was filled with the directions of all 
the events actually detected by the LAT during the six- 
month period. Then, the estimated CR-background map 
was subtracted from that all-particle six-month map to 
produce a residual map that was assumed to contain only 
gamma rays. The number of gamma rays detected dur- 
ing the GRB observation was then calculated by scaling 
the residual map with the ratio of the exposures of the 
two periods (six months exposure over GRB-observation 
exposure). 

The procedures of estimating the CR and gamma-ray 
components of the background described above were re- 
peated for 40 different energy ranges spaced logarith- 
mically from 20 MeV to 300 GeV. For each energy 
range, the CR-background estimations used a different 
dependence of the CR-background rate on the geomag- 
netic coordinates, and 9 and cf) distributions. Similarly, 
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the gamma-ray background estimations used different 
instrument-response functions for the exposure calcu- 
lations. The resulting CR and gamma-ray maps were 
then added to produce all-particle estimated-background 
maps (one for each energy range). The new maps were 
then integrated over the energy-dependent regions of in- 
terest to produce the final all-particle background esti- 
mates. The all-sky fraction of gamma-rays in the back- 
ground appears consistent with the \B\ > 70° gamma- 
ray fraction, such that it contributes ~ 10% towards 
the overall background rate. The results of this method 
were tested against actual LAT data for a variety of du- 
rations, locations in the celestial sphere, energies, and 
data classes. All the distributions of the ratios of the 
estimated over the actually-detected signals followed a 
Gaussian distribution, with width about 15% and cen- 
ter zero (no systematic over- or under-estimation of the 
background) . This accuracy did not depend on the direc- 
tion of the region of interest (e.g. its distance from the 
Galactic plane) and did not have a strong dependence 
(~5%) on the observation's duration. 

2.2.3. Significance calculation 

Using the 'transient' event selection in an energy- 
dependent analysis region centered around the GBM best 
localization and no energy selection, we find 15 events 
(A on ) between To and To + 35 s (time interval where we 
find significant emission in the Nal detectors). From the 
background estimate (Section 2.2.2), the expected num- 
ber of counts in the same region and time interval is 
B C st — 1-3. The gamma-ray contribution to the back- 
ground estimate for this particular burst was about 8%. 
In order to assess carefully the significance of this obser- 
vation, we have used 4 independent statistical methods 
which are described below. 

The first method uses an unbinned likelihood analysis 
of the LAT data which takes into account the energy- 
dependent PSF in an event-by-event basis. This method 
finds a significance of 6.5 comparable to the significances 
found with more simple counting methods described be- 
low. This is due to the fact that our selected events have 
fairly large PSF and thus spatial information within our 
ROI is not that constraining. In the case of GRBs where 
high energy events have been detected, this method is 
expected to provide the highest significance since it fully 
takes into account the spatial information for each event. 

The second method computes the probability of the 
null hypothesis being true (the probability that the ob- 
served number of counts in the on-source region is due to 
a background fluctuation) in a frequentist approach that 
treat s the background u ncertainty in a semi-Bayesian 
way (jConrad et al.l I2003D . Given a certain estimated 
number of background counts B during the on-source 
interval, we compute what is the probability of the ac- 
tual on-source measurement N on being consistent with 
this value. P S up(N on , B) expresses the probability of the 
on-source measurement being equal or superior to N on 
when only statistical fluctuations are considered: 
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estimation method, each possible value for the number 
of background counts is weighted using a Gaussian dis- 
tribution with a mean of B est and a standard deviation 
of 0.15 x B cst (since we estimate our systematics to be 
about 15%): Gauss(B). We then integrate over all the 
possible values of B to compute a weighted probability: 



J^dB Gauss(fl) x P sup (A on ,g) 
r o °° dB Gauss(Ti) 



(2) 



Because of systematic uncertainties in the background 



We applied this method for the numbers mentioned 
above: A on = 15 events, B cs t = 1.3, and we obtained a 
null-hypothesis probability of 1.3 x 10~ 10 , which corre- 
sponds to a significance of 6.6 a. 

The third method is fully Bayesian. The question of 
whether a GRB is detected by the LAT is analyzed as 
an on-source/off-source observation in the time domain. 
For the on-source observation the same parameters are 
used as in the classical (frequentist) analysis, T on = 33.0 
s of live time (corresponding to a clock time of 35.5 s) 
and N on =15. Because the inclination angle of the best 
source position with respect to the instrument boresight 
moved slowly in a 565 s interval around the burst trig- 
ger (with T ff = 525.0 s of live time), we could provide 
an estimate of the off-source background rate. During 
this interval there were N a g = 19 counts detected, and 
the corresponding rate is consistent with the previous 
background estimate B cs t- The spacecraft motion was 
favorable for GRB 080825C, allowing an unusually long 
off-source interval to determine the background rate. For 
some other GRBs, spacecraft motion may cause the back- 
ground in the LAT to vary more quickly, limiting the time 
for which an off-source interval measures the same back- 
ground as that of the on-source interval, and therefore 
limiting the applicability of this method. 

The Bayesian method assumes that the counts during 
the background interval are due to a Poissonian process 
with rate b. To evaluate the probability that a source has 
been detected during the on-source interval, the method 
compares two hypotheses for the Poisson rate during the 
on-source interval: that observed counts are due to the 
same background rate b, or that they are due to a back- 
ground plus source rate b + s. This second hypothesis is 
insufficiently specified to quantitatively solve the prob- 
lem: one must specify some plausible range for the source 
rate s using prior data (i.e., not using the LAT observa- 
tions). We produce a reasonable estimate for the LAT 
source counts by extrapolating the time-integrated spec- 
tral fit of the GBM data to the LAT energy range. The 
photon model from the GBM fit is propagated through 
the LAT response to predict the number of LAT counts. 
A real LAT observation could have a smaller value than 
this estimate because of a statistical fluctuation, or be- 
cause of a spectral break between the GBM and LAT en- 
ergy range. Alternatively, a real LAT observation could 
actually exceed this estimate for the maximum rate if 
there were an additional and distinct spectral component 
to the one found with the GBM data alone. Neverthe- 
less, this is a reasonable 'prior' estimate for the maximum 
counts expected in the LAT, which for GRB 082525C 
is ~60 counts. Un der the assumpti ons described, both 
lLoredo I (p.992) and I Gregory] (|2005[ ) give analytic solu- 
tions for the probability P that the source is detected. 
Using the observational parameters listed above, we find 
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1 — P = 1.2 x 10 -8 which corresponds to a 5.6 a signifi- 
cance. Moreover, this method can provide the probabili- 
ties for the number of events actually originating from the 
source: all 15 events (30%), 14 events (35%), 13 events 
(22%), 12 events (9%), 11 events (3%), <10 events (1%). 

Finally we have computed the significance with a fully 
frequentist method using the on-sour ce/off-source ap - 
proach as described by (equation (17) of lLi &; Ma~lll983[ ). 
This method yields a significance of 6.4 a for the detect- 
ing this burst. 

It should be noted that because such search for LAT 
excess is performed on all GRBs triggered by the GBM 
and other instruments (when the burst is in the LAT 
FoV), it is important to consider multi-trials in our anal- 
ysis. For independent searches as is the case here, the 
post-trials probability threshold for obtaining a 5cr result 

is Ppost-trial = 1 - (1 - P^a) 1/N where N is the number 
of trials and P 5a the 5a probability threshold for a single 
search (~ 5.7 x 1CT 7 ). In the case of GRB 080825C, we 
searched for LAT excess in ~50 bursts triggered by the 
GBM which corresponds to a post-trial probability for 
a 5cr results of P^a, post-trial ~ 1.15 x 10 -8 . This cor- 
responds to a significance of 5.7 which is therefore our 
threshold for a 5<r detection. 

The four independent significance computations pre- 
sented above all yield consistent results similar or above 
this threshold. GRB 080825C is therefore the first G RB 
detected by the LAT instrument (|Bouvier et al.ll2008T ) at 
a high significance level. 

2.2.4. Localization 

Due to th e strong variation of the LAT PSF with pho- 
ton energy (|Atwood et a l. 2009), the on-ground localiza- 
tion of a source depends strongly on its spectral shape. 
For example, 10 photons detected at 100 MeV will yield 
an accuracy of ~1°, while one single photon with an en- 
ergy of 10 GeV will increase the accuracy to ~0.1°. 

The on-ground localization procedure makes use of the 
'transient' class events. It is restricted to the events de- 
tected above 100 MeV (which have a good PSF) in a 15° 
region around the GBM trigger position (RA = 232.2, 
Dec = —4.9). The method is base d on a likeliho o d rati o 
test, following the same steps as iMattox et al.l (|1996j ). 
While these authors used a binned likelihood for the 
analysis of EGRET data, our likelihood function is un- 
binned and uses the instrument PSF on an event basis. It 
also takes into account the various residual backgrounds. 
First, the position of the source and its spectrum are 
left free in the fit, assuming a power-law shape. Then, 
in order to compute an accurate localization error, the 
Test-Statistics (TS, see section [3.1.2l for a complete defi- 
nition) of the point source spectral fit is computed at each 
node of a fine map (5° x 5° for th i s case , with a bin size 
of 0.1°). Following IMattox et all (fl996l) . the TS values 
are interpreted in terms of the chi-squared distribution 
with two degrees of freedom (the two map coordinates). 
Figure [2] shows the obtained error contours around the 
fitted position RA = 233.9, Dec = -4.5, with a 68%, 
90% and 99% statistical error radius of 0.8°, 1.3° and 
2.0°, respectively. 

At large inclination angles, the systematic error of a 
source localization with the LAT is dominated by the 
slight bias in direction reconstruction of low-energy pho- 
tons. This bias is mainly caused by a trigger effect, which 



selects those events that scatter downwards and interact 
with at least three trac ker planes, as requir ed by the in- 
strument trigger logic (|Atwood et al.ll2009h . The bias is 
amplified by the reconstruction efficiency, which is larger 
for tracks near normal incidence. Since the PSF is as- 
sumed to be centered on average on the true directions 
of the gamma rays, this effect translates to a bias in the 
fitted position towards smaller inclination angle. 

This systematic error has been evaluated in two steps, 
using both data and Monte-Carlo simulations. In the 
first step, we studied the LAT performance in localiz- 
ing the Vela pulsar, which is the brightest source and 
has a well-determined position from observations at other 
wavelengths. We found that the fitted position obtained 
from observations where this source was seen at large in- 
clination angles, is biased towards smaller angles. This 
bias is noticeable only when low-energy events (below 
1 GeV) are used and it disappears when high-energy 
events (above 1 GeV) are included in the analysis. The 
agreement found with the prediction of the Monte-Carlo 
simulation allows us to evaluate the bias for any burst 
observation. For instance, Figure [3] shows how the bias 
varies with energy range for a bright gamma-ray point 
source with the same spectral index and inclination an- 
gle as GRB 080916C (|Abdo et al.ll2009al) . 

In the second step, we used the Monte-Carlo simu- 
lation to evaluate the localization systematic error for 
GRB 080825C. We produced a simulation of a point 
source with the same spectral index (~— 2.3) and in- 
clination angle as GRB 080825C (-60°). The position 
obtained with events between 100 MeV and 570 MeV 
showed a deviation of ^0.6° with respect to the input 
position. This systematic error is larger than the one de- 
rived for the very bright GRB 080916C (<0.1°) because 
GRB 080825C occurred at a larger inclination angle and 
had a maximum photon energy well below 1 GeV. 

3. DATA ANALYSIS 
3.1. Temporal analysis 
3.1.1. GBM and LAT light curves 

The Nal, BGO and LAT ('S3' selection) light curves 
are shown in Figure [1] No LAT event is seen in coinci- 
dence with the first bright GBM peak. The first 3 LAT 
events above 80 MeV are detected in a very short time 
window, a few seconds after the GBM trigger, in coinci- 
dence with the second GBM peak. After a quiet period 
where no LAT events are detected up to ~16 s, 4 more 
events are detected within the GBM T90, and another 4 
events after Tgo when the Nal emission has faded close to 
background level. Interestingly the highest energy event, 
with an energy of 572±58 MeV, is detected at ~T +28 s. 

To quantify the different features of this GRB (LAT 
delay, gap, and extended emission), we have performed 
Monte-Carlo simulations of the LAT light curve and es- 
timated the fraction of those simulations that reproduce 
these features. The LAT event distribution was produced 
using Poisson statistics for a constant background rate of 
0.037 Hz (see section |2.2.2[) . an estimated detected sig- 
nal above 80 MeV of 11.7 events (13 events minus 1.3 
estimated background events) and a temporal probabil- 
ity distribution based on the Nal light curve from To to 
To + 35 s. The probability of the different features were 
computed as follow: 
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• delay: the fraction of simulated light curves where 
the first event arrives later than the first actual 
observed photon (at To + 3.252 s). 

• gap: the fraction of simulated light curves which 
include a gap in the middle of the light curve with 
a width larger than the 12.38 s observed. 

• extended emission: the fraction of simulated light 
curves where the last 5 events are detected after 
the timing of the ninth observed event (detected at 
T + 26.570 s). 

This analysis finds weak evidence for the possible de- 
lay or gap features (with chance probabilities of 3.4% and 
0.89%, respectively), but the evidence for temporally ex- 
tended emission in the LAT is more significant, at a 3.7 a 
level. 

3.1.2. High-energy afterglow search 

We searched for possible afterglow emission up to 13 ks 
after the trigger time, but did not find any significant 
emission. Note that the LAT event detected around To + 
47 s is consistent with the expected background event 
every —30 s in the region of interest. 

For the search, we selected the time intervals in which 
the 10° region of interest centered on the LAT location 
(see section 12.2. 4|) was in the LAT field-of-view. The 
burst location exited the field-of-view 1500 s after the 
trigger time, re-entered it — 1 hour later (at Tq + 5200 s) 
for -1600 s, and again -1.5 hour later (at T + 10800 s) 
for -2400 s. 

These periods were split into five time bins with 
roughly exponentially increasing durations : 35.5 s to 
100 s after the trigger time (just at the end of the prompt 
emission), 100 s to 350 s, 350 s to 1500 s, 5200 s to 6800 s 
and finally 10800 s to 13200 s after the trigger time. W e 
used the 'diffuse' event selection (|Atwood et al.l |2009|) . 
which is adapted to faint source studies. 

In each time bin we computed the Test-Statistics (TS) 
of a point source with a power-law spectrum located at 
the position of the burst. Two spectral fits were per- 
formed using the unbinned maximum likelihood method, 
one including only background components (the null hy- 
pothesis), the other also including a test point source 
(the alternative hypothesis), with two possible spectral 
indices: —2.0 (close to the spectral index of the prompt 
emission in the last time bin) and —1.5 (i.e., a harder 
spectrum which could arise from an additional spec- 
tral component at later times). The Galactic diffuse 
emission, the isotro pic diffuse emission (as described in 
lAbdo et al . (2009b)) and two likely blazars close to the 
burst l ocation (OFGL J15 11. 2-0536 and 0FGL J1512.7- 
0905) (|Abdo et all l2009cf ) were included in the back- 
ground model, and their contribution was estimated us- 
ing pre-burst data for the region of interest. The TS, 
which is defined as two times the difference of the log- 
likelihoods between the alternative hypothesis and the 
null hypothesis, was derived in each time bin and for 
each considered spectral index. All TS values are very 
close to zero, indicating a null detection. 

Since no significant afterglow emission was found, 95% 
C.L. upper limits on the flux of the possible emission 
were derived, using a Bayesian method with flat prior 



(iHelene 1 1198"! . This meth od was preferred to the likeli- 
hood profile (|Rolke at ail 120051 ) because of the very low 
count regime. Since the TS is close to zero, the profile is 
not a symmetric parabola and the upper limits derived 
by the likelihood profile method do not have proper sta- 
tistical coverage. The results obtained with the Bayesian 
method are shown in Figure H] and in Table [1] for both 
assumed spectral indices. The upper limit found in the 
first bin is of the same order of magnitude as the flux 
derived at the end of the prompt emission. Then the up- 
per limit decreases when the considered observation time 
increases. 

The same method was used to derive 95% C.L. upper 
limits on the flux in the prompt emission phase where 
the LAT did not detect any photon (time bins (a) and 
(c)). We assumed a spectral index extrapolating the 
GBM spectrum. These results are shown in Figure [4] 
and Tabled] as well. 

Between bin (e) and bin (h) we assumed that the flux 
decreases with time like t a . We used the fitted flux values 
on each time bin, and for bins with TS zero, we set the 
flux to an arbitrary low value (this choice has no impact 
on the following result). The 68% C.L. error bars were 
determined by the likelihood profile method. Although 
its coverage is not correct for the bins of TS zero, this 
error bar is accurate enough for the present purpose of fit- 
ting the flux decay. It also ensures homogeneity among 
the error bars used for the fit. Since the fitted decay 
slope was not sign ificant, we used a \ 2 profile method 
(Rol ke at al.l l2005) to derive upper limits on the flux de- 
cay slope for both assumed spectral indices : a < —2.08 
(95% C.L.) for a spectral index of —1.5 and a < —1.77 
(95% C.L.) for a spectral index of —2.0. 

3.2. Time-resolved spectroscopy 

We have performed detailed spectroscopy of the com- 
bined GBM and LAT data, for the whole duration of 
the burst and also time-resolved analysis based on the 
temporal structures observed in the GBM and LAT light 
curves. Figure [T] shows the five time intervals (a to e) 
adopted for this analysis. The fits were performed with 
the spectral analysis software package RMFIT (version 
2.5), using log- likelihood as the fitting statistic since chi- 
squared is less appropriate due to the low number of 
events at high gamma-ray energies in this burst. The 
variable GBM background is subtracted for all detectors 
individually by fitting an energy-dependent, second or- 
der polynomial to background data of —300 s before and 
-300 s after the GRB. For t he LA T data, the 'S3' event 
selection is used (see section |2".2. 1[) . We used a LAT en- 
ergy range from 80 MeV to 200 GeV (since energy bins 
with no detection still contain useful information) and 
adopted a constant but energy-dependent background 
rate, as described in section [2. 2. 21 For the GBM data we 
used the standard 128 energy bins of the CSPEC data- 
type, but only the channels above 8 keV in the Nal detec- 
tors, above 240 keV in the BGO detectors, and rejecting 
the overflow channels in both Nals and BGOs. 

We considered different empirical models in the spec- 
tral analysis: a simple power law, a power law with 
a high-energy ex ponential cutoff, a Band function 
(|Band et al.l Il993f ) which smoothly connects two power 
laws, a Band function with a high-energy exponential 
cutoff, and a Band function with an additional power law. 
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The high-energy exponential cutoff was implemented by 
multiplying the original spectrum by exp(— E/E cuto s). 
The main results of our combined GBM and LAT anal- 
ysis are shown in Tables [2] (spectral parameters) and [3] 
(flux and fluence in 50 - 300 keV and 100 - 600 MeV), 
for which we used the LAT data and responses obtained 
with the 'S3' cut. The time- integrated spectrum (shown 
in Figure [5]) and the time bins (b) to (d) are best fit 
to a significant degree by the Band function, and thus 
only the best fit parameters and associated statistical 
uncertainties of the Band function are provided in Ta- 
ble [5] for these time bins. The best fit parameters for the 
Band function are also given for time bin (a), although 
an addition of a possible exponential cutoff is discussed 
later in this section. For time interval (e), however, the 
spectrum is adequately described by a single power law, 
and adding more parameters does not improve the fit. 
Adding a power law function to the Band function does 
not improve the time-integrated and time-resolved fits. 
Besides spectral fits to the whole GBM & LAT energy 
range, we also fitted the GBM data alone and propagated 
the photon model from the GBM fit through the LAT re- 
sponse to obtain the predicted numbers of LAT counts. 
A comparison between the expected and observed num- 
bers of events shows that there is no need for an extra 
emission component besides the Band function (or power 
law for time bin (e)). 

Time bins (a) to (d) display the ty pical hard-to-soft 
evolution of -E pC ak (N orris et al J 1 198a ). which is the en- 
ergy at which the Band function peaks in vF v spectrum, 
starting at almost 300 keV and decreasing to ~150 keV. 
Except for the second interval, the values of the low- 
energy spectral index, a, and the high-energy spectral 
index, j3, are constant within their uncertainties. The 
evolution of the spectral parameters, flux, fluence, and 
the flux ratio between the two energy ranges, are shown 
in Figure O The best-fit model spectra for time bins (a) 
to (e) are shown in Figure [7] along with their 1 a confi- 
dence intervals. 

A significant hardening of the spectrum at high energy 
is observed after ~25 sec. From the light curves (Fig- 
ure []} it is clear that while the GBM emission in time bin 
(e) is just above the background level, there is significant 
emission in the LAT, including the two highest energy 
events detected for this GRB. This spectral hardening is 
also reflected in the power-law index of —1.95 ±0.05 for 
time bin (e) (see Table [2J , which is significantly harder 
than the values of (3 in the earlier time bins. Looking 
at the evolution of Speak, one could argue that the spec- 
trum is affected by curvature at the low-energy end and 
that a Band function is not preferred over a single power 
law due to poor statistics. However, a spectral fit in in- 
terval (e) of all the data above 300 keV gives a softer but 
consistent power-law index of —2.10 ± 0.08. This (3- value 
is closer to the j3- values found for interval (a) to (d), but 
still significantly harder. 

We have searched for possible departures from a simple 
Band function in the different time bins by performing 
a likelihood ratio test comparing a simple Band func- 
tion with a Band function multiplied by an exponential 
cutoff term, exp(— E/E cuto e)- A significance of 4.3 er was 
found for an exponential cutoff in time bin (a) with a cut- 
off energy around E cuto s = l-77^o.56 MeV (with follow- 



ing Band function parameters: a ~ —0.57, /3 ~ —1.64, 
Epeak ~ 211 keV). We investigated the dependence of 
this significance with the systematics of our instruments 
and found the strongest effect to be a ± 15% variation in 
the BGO effective area which can bring the significance 
down to ~3.7<7. With 5 time bins, this is not strong 
enough to claim the existence of an exponential cutoff. 

We have performed a careful investigation of the effect 
of systematic uncertainties on the spectral parameters 
of GRB 080825C. We considered the following systemat- 
ics in the LAT, BGO and Nal detectors: effective area, 
energy dispersion and background subtraction. For the 
LAT the methodology we use is to modify the Instru- 
ment Response Functions (IRFs) in order to take into 
account those various effects. Systematics on the LAT 
effective area have been derived from a study of the Vela 
pulsar and Earth albedo photons because it is possible 
in both cases to extract an extremely pure gamma-ray 
sample. Using this sample, the uncertainty on the LAT 
effective area has been computed as a function of en- 
ergy to be 10% below 100 MeV, 5% around 1 GeV and 
20% above 10 GeV. Based on this estimate, some spe- 
cial IRFs have been created to encompass the extreme 
scenarios in overall normalization and slope of the LAT 
effective area. The uncertainty in energy measurement 
for the LAT is estimated to be of the order of 5%. We 
adopted a 10% uncertainty in the Nal and BGO effec- 
tive area (both overall normalization and slope). Finally 
we also considered the uncertainty in the choice of off- 
timing sample for Nal and BGO background subtraction 
and found the corresponding systematics to be negligi- 
ble. In each time bin, the error values on the spectral 
parameters are found to be similar to or smaller than 
the statistical uncertainties reported in Table [2j except 
for the case of time bin (a) and (b) where the systemat- 
ics on -Epoak and the normalization are found to be about 
two times and three times larger, respectively. Table 
reports the systematics found for each parameter of the 
time-integrated best fit Band function reported in Ta- 
ble [2] as well as the predominant systematic effects. We 
note that similar trends for t he systematic uncer tainties 
were found for GRB 0809 16C IjAbdo et alJl2009al) . which 
had many more LAT events than GRB 080825C. 

4. DISCUSSION 

4.1. Theoretical interpretation of GRB 080825 C 
observations 

GRB 080825C is the first GRB detected by the Fermi 
LAT. Furthermore, it is the first GRB to show a hint 
of a time delay between the onset of the high-energy (> 
100 MeV) emission relative to the low-energy (sub-MeV) 
emission (later see n much more clearly in GRB 080916C; 
lAbdo et aLll2009ah . Moreover, there appears to be a lo- 
cal minimum in the high-energy flux between the initial 
narrow (< 1 s) LAT spike, at —3 - 4 s after the GBM 
trigger, and the later, much broader peak (between ~16 s 
and ~31 s). GRB 080825C is the first GRB to show this 
possible feature in the high-energy emission. 

The late broad peak in the high-energy emission has a 
duration comparable to its peak time as measured from 
the GRB trigger time, and these two timescales are also 
similar to the duration of the low-energy emission (which 
has Tg ~ 27 s and T 50 ~ 13 s). The fact that these three 
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timescales are comparable is naturally accounted for if 
this late broad peak signifies the onset of the emission 
from the external (reverse or forward) shocks, if the re- 
verse shoc k is (at least mildly) relativistic (see Figure 1 of 
ISarllT997h . This occurs around the deceleration time, as 
the GRB outflow is decelerated by a reverse shock while 
it drives a highly-relativistic forward shock into the ex- 
ternal medium. In this scenario the late broad LAT peak 
can be either synchrotron self-Compton (SSC) emission 
from the reverse shock or external-Compton (EC) emis- 
sion from the forward-reverse shock system that forms as 
the GRB outflow is decelerated by th e external medium 
(as was suggested for GRB 941017; iGranot fe Guettal 
l200llPe'er fc Waxman|[200l . In the latter case the seed 
synchrotron photons can come from the reverse shock 
and the upscattering relat ivistic electrons fr om the for- 
ward shock, or vice versa (|Wang et al.H200lD . 

In time bin (e) , which contains the second half of the 
broad high-energy LAT peak while the GBM emission is 
almost back to background, the spectrum is well fit by 
a single power law with a photon index of —1.95 ± 0.05, 
which is significantly harder than the values of the high- 
energy photon index (3 in all of the earlier time bins. 
Moreover, the two highest energy photons in this GRB 
are in time bin (e). This may suggest that the late time 
wide high-energy peak arises from a separate spectral 
component to that responsible for the low-energy emis- 
sion, which is consistent with an origin from a distinct 
physical region (and in particular with an external shock 
origin, as mentioned above). If there is still some con- 
tribution from the low-energy spectral component, then 
the photon index of the high-energy spectral component, 
which dominates in time bin (e) but is sub-dominant ear- 
lier on, could be somewhat harder than the measured 
value. A rather flat vF v spectrum around ^100 MeV 
might suggest that the peak photon energy of this spec- 
tral component is close to this energy range, which may 
be readily obtained for SSC from the reverse shock, and 
possibly for the EC process described above. Such a 
spectrum would disfavor SSC from the forward shock, 
which typicall y peaks at TeV energies near the decel- 
eration time (I Wang et al.ll27M lG7an"ot fc Guettal 120031 : 
IPe'er fc Waxmanll2004n . 

The rate of the rise in flux up to the wide peak is 
hard to quantify, and therefore does not significantly con- 
strain the theoretical models. The upper limits on the 
late time high-energy flux imply a flux decay at least as 
steep as ~ i -1 . This is consistent with the steep decay 
that is expected from the reverse shock emission after 
the passage of the reverse shock, but only b arely consis- 
tent with forward shock SSC emission (e.g.. iSari fc Esir] 
I2001H for fast cooling for a sufficiently soft electron en- 
ergy distribution, corresponding to a power-law index of 
p > 2.6 — 2.7. Other emission mechanisms might also be 
able to produce a sufficiently steep flux decay rate. 

The possible high-energy cutoff around E cuto e ~ 
1.8 MeV in time bin (a), i.e. the first 2.7 s from 
the trigger time, if true, is very interesting, as it does 
not appear to have a good simple explanation. The 
lack of a redshift measurement for GRB 080825C, how- 
ever, makes it difficult to draw very strong conclusions, 
even in this case. In time bin (a), -E pea k ~ 290 keV 
is a factor of ~6 lower than -E cu toff, thus disfavoring 
a quasi-thermal spectrum. A possible origin of such 



a cutoff from in trinsic opacity to pair product i on in 
the source Ce.g.. iLithwick fc Sari I [2 001: Ba riTjel 120061 
IGranot. Cohen- Tanugi fc do Couto e Silval [20081 ) might 
be hard to reconcile with the non-thermal spectrum and 
sharp peak in the GBM light curve in time bin (a), which 
both imply that the Thomson optical depth tt of the 
pairs produced in the source cannot be 3> 1, for any rea- 
sonable redshift (that is not C 1). Such a reasonable red- 
shift would, in turn, typically imply a bulk Lorentz factor 
of the emitting region of T > 100, and make it very dif- 
ficult to produce -E C utoff ~ 1.8 MeV, corresponding to a 
comoving energy of E' mtoB ~ 18(1 + z)(r/100) _1 keV < 
m e c 2 . This can be understood as follows. Producing the 
observed cutoff through intrinsic opacity to pair produc- 
tion requires T-y^(E' cutoS ) ~ 1. Since the opacity to pair- 
production for a photon of energy E' cutoS * s produced 
mainly by photons of energy ~ E' an ~ (m e c 2 ) 2 / E' cntoS 
(i.e. near threshold), and the cross-section to pair pro- 
duction near threshold is of the order of the Thomson 
cross-section (ot), then r 77 (£") ~ TT[(m e c 2 ) 2 / E'] where 
tt{E') is the Thomson optical depth of the pairs that are 
produced if all the photons of energy near or above E 1 
pair produce. In particular, taking E' — E' cutoS implies 
that TT{E' an ) ~ 1 (i.e. the Thomson optical depth of the 
pairs produced by all the photons of energy E' > E' an is 
of order unity). In our scenario r 77 (E' cutoS ) ~ 1 implies 
that all photons above E' cutoS pair produce, and since 
they are much more numerous than the photons above 
E' an this produces a very large Thomson optical depth 
in pairs, r T (E' cutoS ) ~ r 77 (^ n ) ~ {E' & J E'^)- 1 ^ ~ 

{E' cutoS /m e c 2 ) 2(1+ ^ > 1 (where < -1 is the high- 
energy photon index). This is inconsistent with the ob- 
served spectrum and light curve, as mentioned above, 
since such a large optical depth would thermalize the 
spectrum and suppress the temporal variability. Alter- 
natively, a high-energy cutoff may reflect the energy spec- 
trum of the accelerated electrons, but it is not clear why 
a power-law would extend over a very narrow range in 
this case and a much larger dynamical range in most 
other cases. 

4.2. Comparison to recent GRBs with high-energy 
emission 

GRB 080514B was detected by AGILE, triggering on 
board both the SuperAGILE X-ray detector (SA) and 
the AGILE/MCAL detector, and was observed by th e 
GRID instrument up to 300 MeV ([Giuliani et al.l f2008). 
The gamma-ray data above 30 MeV show a significant 
extended duration with respect to emission in the hard 
X-ray and soft gamma-ray energy bands. A detailed 
analysis of this component was not possible due to lim- 
ited statistics. Another interesting fact is that in the case 
of GRB 080514B, some GRID events are detected when 
the low-energy emission has faded beyond detectability 
both in MCAL and SA instruments. Similarly, in the 
case of GRB 080825C, emission in the Nal detectors has 
receded to background level when the last LAT events are 
detected. In both cases, the behavior of the high-energy 
emission in the last part of the prompt emission does 
not seem to be correlated with that of the low-energy 
emission. This suggests that it may originate from a 
physically distinct emission region. The similar behavior 
observed in the high-energy emission of GRBs 080514B 
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and 080825C may suggest a common explanation. 

The longer duration of the high-energy (> 100 MeV) 
emission relative to the low-energy (< 1 MeV) emis- 
sion, and to some extent (though with much smaller 
statistical significance) its later onset, are present not 
only in GRB 080825C, but al so in the two subse quent 
LAT GRBs, GRB 080916C (lAbdo et all l2009ah and 
GRB 081024B (|Omodei et al.ll2008[) . Therefore, they ap- 
pear to be common features of at least the first three LAT 
GRBs. Additionally, the onset of the LAT high-energy 
emission coincides with the second peak in the GBM low- 
energy li ght curve in all the se three LAT GRBs . This 
has led lAbdo et al.l (|2009a) to favor an interpretation 
in which the first and second spikes of the low-energy 
light curve originate from two distinct emission regions 
with different intrinsic emission spectra, where the sec- 
ond spike has a harder high-energy photon index f3. This 
therefore accounts for the LAT detection of the second 
spike and non-detection of the first spike. This interpre- 
tation is further supported in the case of GRB 0809 16C 
by the clear change in both the low-energy (a) and high- 
energy (/3) photon indices between the first and second 
peaks of the GBM (low-energy) light curve. This is 
not clearly seen in GRB 080825C, possibly since it was 
not nearly as bright as GRB 0809 16C. Other explana- 
tions for the delayed onset are p ossible, and s e veral o ther 
models were also considered by Ab do et al.l (|2009a[ ) for 
GRB 080916C. It was hard to clearly distinguish between 
the different models on the basis of the available data, 
despite the extreme brightness of GRB 080916C (which 
had an extremely large fluence of 2.4 x 10 -4 erg cm~ 2 
and a record breaking E lMo of 8.8 x 10 54 erg). For 
GRB 080825C it is even harder to distinguish between 
the different possible explanations. 

Temporally extended high-energy (> 100 MeV) 
emission, which lasted longer than the prompt low- 
energy (sub-MeV) emission, was detected not only 
in GRB 080825C, but also in the other two LAT 
GRBs mentioned above (080916C and 081024B), as 
well as in the AGI LE GRB q80514 B and the 
EGRET GRB 940217 (jHurlev et al.|[i~99l . However, 
GRB 080825C is the only one so far that shows a hint of 
a minimum in the high-energy flux between the early 
and late high-energy emission, which strengthens the 
case for an origin from a distinct physical region, as 
discussed above. In GRB 080825C the late time high- 
energy emission has a harder photon index than the ear- 
lier high-energy emission, which is consistent with an 
external shock origin. The opposite is true for GRB 
080916C, which favors a somewhat different model such 
as external Compton scattering of late time X-ray flare 
photons by forward shock electrons. In GRB 080825C 
the highest-energy photon was detected after 28 s when 
the low-energy emission was already down to a very low 
level, while in GRB 080916C the highest energy pho- 
ton was detected after 17 s, while the low-energy emis- 
sion was still very active. Moreover, the observed dura- 
tion of the high-energy emission in GRB 080825C lasted 
only slightly longer than the low energy emission (up to 
~31 s compared to T 90 ~ 27 s) while in GRB 080916C 
it lasted much longer (more than ~10 3 s, compared 
to T90 = 66 s). While the high-energy emission in 
GRB 080916C was detected for a longer time, partly be- 



cause it was generally much brighter than GRB 080825C, 
it is also possibly because the flux decayed more slowly at 
late times (as i _1 - 2±0 - 2 , compared to steeper than ~ f -1 - 7 
in GRB 080825C). 

5. CONCLUSIONS 

GRB 080825C is the first GRB detected by the LAT, 
with 13 events above 80 MeV and a detection significance 
of ^6 a. The highest energy events, up to ~600 MeV, are 
detected at late times, ~25 — 35 s after the GBM trigger, 
when the emission in GBM has decreased close to the 
background level. The lack of > 1 GeV events in the LAT 
and the large angle of the source to the LAT boresight 
result in a localization uncertainty of 1.1° (statistical plus 
systematic) at the 1 a level. 

The prompt emission spectrum from both instruments 
onboard Fermi covers over 5 decades in energy. We have 
performed time-resolved spectral analysis using the two 
GBM Nal detectors with the brightest GRB signal, both 
GBM BGO detectors, and for the LAT with an event 
selection scheme that is optimized for GRB analysis. We 
have carefully taken the energy-dependent backgrounds 
into account for both GBM and LAT, and studied the 
systematic uncertainties in the spectral analysis. The 
time-integrated and time-resolved spectra are well fit by 
the Band function with a hard-to-soft evolution in the 
first 25 s: E vea y evolves from ^300 to ~150 keV, the 
high-energy power-law index /? is constant at a value of 
~— 2.5, while the low-energy power-law index a is fairly 
constant except for the second time bin which contains 
the first LAT events. In the last time bin, ~25 — 35 s after 
the GBM trigger time, the GBM data are barely above 
background level, and the spectrum is best fit by a single 
power law with an index of ~— 2 which is significantly 
harder than the (3 values of the earlier intervals. 

The duration and start time of the late broad peak 
in the high-energy emission, ~16 — 31 s after the trig- 
ger, suggest that this peak is emitted by the external 
reverse or forward shocks, rather than by internal dissi- 
pation within the GRB outflow (e.g., internal shocks or 
magnetic reconnection) . The relatively fast flux decay 
after this peak slightly favors a reverse-forward shock 
'external' Compton origin over a forward shock SSC ori- 
gin. Although the origin and emission mechanism for this 
late peak cannot be conclusively determined because of 
low number statistics (and the lack of observations at 
X-ray or optical wavelengths, due to the poor GRB lo- 
calization) , the external shock origin is further supported 
by the change in spectral behavior, in particular of the 
spectral index, at these late times. Observations of more, 
brighter GRBs with both GBM and LAT will be able to 
test this hypothesis. 
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TABLE 1 

95% C.L. UPPER LIMIT ON THE PHOTON FLUX IN THE LAT ENERGY RANGE DERIVED BY A BAYESIAN METHOD WITH A FLAT PRIOR. In TIME 
BINS (A) AND (C) OF THE PROMPT EMISSION THE CONTINUATION OF THE SPECTRUM FITTED ON GBM DATA WAS ASSUMED. In TIME BINS (f) 

to (j) two different cases were studied: a spectral index of —2.0 (close to the spectral index of the prompt emission in 
the last time bin) and of —1.5 (i.e. a harder spectrum which could arise from an additional spectral component at later 

times) . 



Time 


Spectral 


Photon flux above 100 MeV 


Bins (s) 


index 


95% C.L. upper limit (10~ 6 ph cm" 2 s" 1 ) 


a : 0.00 - 2.69 


-2.54 


1750 


c : 4.74 - 12.93 


-2.62 


610 


f : 35.5 - 100. 


-2.0 / -1.5 


170 / 110 


g : 100. - 350. 


-2.0 / -1.5 


14 / 9.1 


h : 350. - 1500. 


-2.0 / -1.5 


2.7 / 1.5 


i : 5200. - 6800. 


-2.0 / -1.5 


1.2 / 0.79 


j : 10800. - 13200. 


-2.0 / -1.5 


1.1 / 0.75 



TABLE 2 

Time-integrated and time-resolved spectral analysis results for GRB 080825C. Band function best fit parameters are 

PROVIDED FOR ALL SPECTRA, EXCEPT FOR TIME INTERVAL (e) WHICH IS ADEQUATELY FIT BY A SINGLE POWER LAW. 



Time Range (s) 


A (10~ a 7 cm" 2 s" 1 keV -1 ) 


a 


P 


£ pcak (keV) 


a: 0.00 - 2.69 
b: 2.69 - 4.74 
c: 4.74 - 12.93 
d: 12.93 - 25.22 
e: 25.22 - 35.46 


75l b 5 

1381" 
44 ±4 
47 ±1 

1.2 ±0.1 (at 100 keV) 


-0.76 ± 0.05 
-0.52 ± 0.06 
-0.81 ± 0.06 
-0 72 + 07 

u - lz -0.06 

N.A. 


-2 54+ lhn 
Z ' 04 -0.17 

r, o 7 +0.06 

-0.08 

-2 62 + 14 
z,uz -0.25 

-2 45+ ' 07 
-0.10 

-1.95 ±0.05 


291 1% 
2101 14 , 
183 ± 13 
152 ±9 
N.A. 


0.00 - 35.46 


37 ±2 


-0.79 ± 0.03 


_2 42+ U ' U4 
Z - 4Z -0.05 


198 ±8 



TABLE 3 

Flux and fluence in the 50 - 300 keV and 100 - 600 MeV energy ranges for the time intervals and spectral parameters 

presented in table [2] the upper limits are given at the 2 a level. 



Time 
bin 


Flux 50-300 keV 
(10~ 7 erg s -1 cm -2 ) 


Flux 100-600 MeV 
(10 -8 erg s _1 cm" 2 ) 


Fluence 50 -300keV 
(10~ 6 erg cm -2 ) 


Fluence 100-6OO MeV 
(10~ 7 erg cm~ 2 ) 


a 


16.2 ± 0.2 


< 20.1 


4.35 ± 0.05 


< 5.39 


b 


20.9 ±0.3 


21.1 ± 7.2 


4.28 ± 0.06 


4.32 ± 1.47 


c 


6.57 ±0.11 


< 3.78 


5.38 ± 0.09 


< 3.09 


d 


5.50 ±0.07 


2.50 ±0.97 


6.76 ± 0.09 


3.07 ± 1.19 


e 


0.34 ±0.03 


5.03 ± 1.68 


0.35 ± 0.03 


5.15 ± 1.72 


Full 


5.98 ±0.04 


4.35 ±0.88 


21.2 ±0.1 


15.4 ±3.1 



TABLE 4 

Systematic uncertainties as well as the predominant effect for each spectral parameter of the time-integrated best fit 

Band function reported in Table [2] 





Norm 


a 


P 


Speak (keV) 


systematic error 


±15% 


±0.03 


±0.03 


±8 


predominant systematics 


Nal eff. area 


Nal off. area 


BGO eff. area 


Nal &c BGO eff. area 
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Time since GBM Trigger (s 



Fig. 1. — Light curves of GRB 080825C observed by the GBM (Nal & BGO) and LAT instruments; top two panels arc background 
subtracted. The LAT light curve has been generated using events which passed the 'S3' event selection above 80 MeV (which are also the 
events used for our spectral analysis). Black dots, along with their error bars (systematic uncertainty in the LAT energy measurement) 
represent the 1 a energy range (right y-axis) for each LAT event. The vertical dash-dotted lines indicate the time bins used in our 
time-resolved spectral analysis. 
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Fig. 2. — LAT on-ground localization for GRB 080825C: RA=233.9, Dec=-4.5. The contours show the containment regions for confidence 
levels: 68%, 90%, and 99%. Equivalent containment radii can be computed: 0.8° (68%), 1.3° (90%) and 2.0° (99%). 
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Fig. 3. — Dependence of the localization bias on the energy range of the events used. These results are obtained from a simulation of 
LAT observations of a point source similar to GRB 080916C (same slope and inclination angle). The use of low-energy events only (below 
1 GeV) yields a noticeable bias, while no bias is seen when using events only above 1 GeV. 
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Fig. 4. — Squares: photon fluxes (or 95 % C.L. upper limits) measured during the prompt emission, in the five time bins defined previously. 
Filled circles (resp. open circles): 95% C.L. upper limits on the afterglow emission photon flux in each considered time bin, assuming a 
power-law spectral shape of index —2.0 (resp. —1.5). 
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Fig. 5.— Time-integrated (T + 0.004 to T + 35.456 s) count spectrum of GRB 080825C of the GBM (Nal and BGO) and LAT data. 
The spectrum is well fit by a Band function spanning ~5 decades of energy. LAT data has been separated into 'FRONT' and 'BAC K' data 
sets which respectively correspond to events converted in the upper and lower part of the tracker instrument (Atwood ct al. 2009). 
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Fig. 6. — Left: Temporal evolution of the spectral parameters for GRB 080825C: the low-energy index a (top), the high-energy index fi 
(middle), and the peak energy i? pca k (bottom). The last time bin is adequately fit by a single power-law function, of which the index is 
plotted in the middle panel. Right: Temporal evolution of the flux (top) and fluence (middle) in two energy ranges: 50 — 300 keV (solid 
circles) and 100 - 600 MeV (open circles); and of the ratio (bottom) between the high-energy and low-energy flux. The upper limits are 
given at the 2 a level. 
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Energy (keV) 

Fig. 7. — The best-fit spectra for time bins (a) to (e) are shown in thick solid lines that reach up to the largest detected photon energy 
in each time bin, while the corresponding (same color) thin dashed lines represent the 1 a confidence contours for each fit. Time bins (a) 
to (d) are fit by a Band function and time bin (e) is fit by a power-law spectral model. 



